Shape resonances in the superconducting order parameter of ultrathin nanowires 
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We study the shape resonance effect associated with the confined transverse superconducting 
modes of a cylindrical nanowire in the clean limit. Results of numerical investigations of the 
Bogoliubov-de Gennes equations show significant deviations of the energy gap parameter from its 
bulk value with a profound effect on the transition temperature. The most striking is that the size 
of the resonances is found to be by about order of magnitude larger than in ultrathin metallic films 
with the same width. 
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^ , Modern rapid miniaturization of electronic circuits requires good understanding of basic mechanisms responsible 
5h ' for the electronic properties of nanoscale structures. The most important point about these structures is that the 
Or quantum-confinement effects play the corner-stone role in this case. One can even say in general that recent success in 
nanofabrication technique has resulted in great interest in various artificial physical systems with unusual phenomena 
driven by the quantum confinement (quantum dots, nanoscale semiconductors, nanosuperconductors, etc.). The 
C$ i quantum-confined superconductivity is here of special interest due to the macroscopic quantum character: any effect 
on electron wave functions manifests itself directly in the superconducting order parameter. 

An obvious consequence of the confinement in a nanoscale superconducting structure is nonuniform spatial dis- 
tribution of the superconducting condensate because, as it is known since the classical papers by Gor'kov Q and 
Bogoliubov 0, the superconducting order parameter can be interpreted as the wave function of the center-of-mass 
motion of a Cooper pair. It is also known that the Cooper-pair wave function involves important in- medium terms 0. 
In the presence of the electron confinement these terms can result in shape resonances in the energy gap parameter, 
another confinement effect first found and investigated in the paper by Blatt and Thompson Q for ultrathin metallic 
films. A shape resonance in the dependence of the energy-gap parameter on the specimen dimensions can occur any 
, time when an electron subband appearing due to the size quantization passes through the Fermi surface Q. Stronj 
indications for such behaviour are found not only in ultrathin films but also in nanoparticles (see, for example, Refs. [5^ 
and and superfluid nuclei 0,0 (as it had been predicted by Blatt and Thompson). Recently a new technique of 
electrodeposition into extended nanopores has been developed 0, which makes it possible to produce single-crystal 
nanowires of high quality. Thus, the shape resonances in the nanowire superconducting order parameter can be 
investigated in the clean limit, with a direct link to the microscopic (BSC) theor y. I n particular, it is of importance 
to explore the situation where the QPS (quantum phase slips) regime |^. Iiot ITU Il2j|) is expected to generate a new 
low-tcmpcrature metallic phase with proliferating quantum phase slips of the superconducting order parameter (for 
radii less than 5 nm) |13|. The point is that when calculating the QPS energy barrier Il.'il| in ultrathin nanowires, 
one assumes the superconducting order parameter being uniform in the transverse direction but with the phase slips 
in the longitudinal one 01 . The absolute value of the order parameter is set to be equal to the bulk one. However, 
at a resonant point the superconducting condensate shows significant spatial variations in the transverse direction, 
and its mean absolute value can be much larger than that of the bulk material. Note that the quantum confinement 
can influence the superconducting state by two channels. The first is due to change in electron wave functions (this 
is a main reason for the shape resonance effect). Whereas the second is connected with the confined phonons. The 
both channels were examined for ultrathin metallic films [T^j . which makes it possible to expect that confinement 
modifications of the phonon modes can produce quantitative corrections if the nanowire width is less than or about 
2 nm. In the present work the first channel is under consideration, while the phonons are taken to be the same as 
in the bulk material. Thus, below we investigate the shape resonance effect associated with the confined transverse 
superconducting modes of an ultrathin nanowire in the clean limit. 

To explore the superconducting order parameter varying with position, one should use the Bogoliubov-de Gennes 
(BdG) equations We are interested in numerical solutions of BdG equations taken in the absence of magnetic 
field for a superconducting cylinder with the radius R and length L. In all the calculations L remains the same and 
equal to 2000 nm while R is varied from 1 nm to 10 nm. In the absence of magnetic field the superconducting order 
parameter A(r) can be chosen as a real quantity (phase effects are beyond the scope of our consideration) and the 
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FIG. 1: The relative gap of the quasiparticle spectrum An/Atuik (panels (al) and (bl)) and relative chemical potential Ha/nbuik 
(panels (a2) and (b2)) versus the nanowire radius R: the panels (al) and (a2) represent the data calculated for n e — 3.878 nm~ 3 
(the bulk chemical potential fibuik = 900 meV); the panel (bl) and (b2) are the data for n e = 20 nm~ 3 (fibuik = 2687 meV). 
Squares are the results of numerical investigations of the BdG equations, the solid line is the spline interpolation. 



BdG equations have the form 

E n u n {r) = (-t^V 2 - Aiju„(r) + A(r)w„(r), (1) 

E n v n (r) = A(r)u n (r) - (-^V 2 - fj)v n (r), (2) 
V 2m* / 

where E n stands for the quasiparticle spectrum, /i is the chemical potential and m* denotes the electron band mass. 
The single electron wave functions u n and v n make a contribution to the order parameter via the self-consistency 
relation 

A(r)= fl 5]tt B (rK(r)(l-2/(^)) ) (3) 

n 

where g is the coupling constant and f(x) is the Fermi function f(x) — l/(exp([3x) + 1), f3 = l/(fcsT) with T the 
temperature and fee the Boltzmann constant. Summation in Eq. © is taken over the eigenstates with the kinetic 
energy (including the chemical potential) within the window [-fiuo,^], and u>d is the Debye frequency. The 
chemical potential is fixed by 

n e = 2j2 [Mr)| 2 /(£„) + K(r)| 2 (l - /(£„))] (4) 
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FIG. 2: The superconducting order parameter A(p) versus p at n e = 3.878 nm and T = 0: (a), the resonance points 
R = 5.06 nm (1) and R = 5.75 nm (2) for the Al parametric set; (b), the resonance point R = 9.642 nm for the Sn parameters. 



with n e the total electron density. Introducing the cylindrical coordinates p, ip, z, we can write A(r) = A(p) as the 
periodical boundary conditions are implied in the longitudinal (z) direction. In this case we get 



u jm k(r) = u jmk (p) —j= —j=, v jmk (r) = v jmk (p) — = —j=, 



(5) 



where n — (j, m, k) with j the quantum number associated with p— coordinate, to the azimuthal quantum number 
and k the wave vector in z-direction. Substituting Eq. J5J into Eqs. Q and J5J , we recast the BdG equations in terms 
of u and v: 



Ejmk Uj m fc (/?) 



/ d 2 i c m° , 2 "\ 
\do 2 odo d 1 / 



1 d 



2m* \dp 2 p dp p 



Ujmk(p) + A(p)v jmk (p), 



Ejmk Vj m k(p) = A(p) u jmk{p) 



h^r + -^ ^-fc 2 - 



1 d 



2m* \ dp 2 p dp p 2 
Due to the electron confinement in the transverse directions we should put 

U Jmk (R) = Vj m k(R) = 0. 



Vjmk(p)- 



(6) 
(7) 

(8) 



Then, u— and v— functions are expanded in terms of the Bcssel functions, and Eqs. © and J7J) can be converted into 
a matrix form convenient for the numerical iteration procedure typical for a self-consistent mean- field treatment. 

Numerical investigation of Eqs. (JSJ and Q reveals the resonance structure of the energy gap parameter of a 
cylindrical nanowire in agreement with the results for ultrathin films and nanoparticles. But what is surprising is that 
the resonance size is now by order of magnitude larger than in ultrathin films of the same width (see the results for 
ultrathin metallic films with the width 10 — 20 A 4]). In Fig.l the R— dependent gap of the quasiparticle spectrum 
An together with the chemical potential pu are given for 0.8 nm < R < 1.8 nm at the zero temperature. The 
Al— parameter set was used in the calculations: gN(Q) = 0.18 (N(0) is the bulk energy density of states for one spin 
projection at the Fermi surface) and Huo/kB — 375 °K 0. The electron band mass was set to the free electron mass. 
To have a feeling about influence of the total electron density on the shape resonances, we performed calculations 
for the densities n e — 3.878 nm -3 (with the bulk chemical potential pbuik = 900 meV) and n e = 20 nm~ 3 (with 
Pbuik = 2867 meV). Any time when the number of electron subbands (associated with the quantum numbers j and m) 
below the Fermi level gets less by 1 or 2 with decrease of the wire radius R (1 for the situation of j dropping, 2 for the 
case of | to | decreasing), a resonance appears in the energy gap dependence on R. The same manifests itself in a sharp 
change of the chemical potential derivative. As it is seen, there is no essential difference in resonance magnitudes of 
panels (al) and (bl) of FigQ] However, the larger electron density, the larger resonance density: the distance between 
two neighboring resonances is roughly proportional to 1/kp with kp the Fermi wavenumber of the bulk material. It is 
interesting that the resonance effect changes the bulk BCS relation A{, It ifc/(fcsT c ) w 1.76 with T c the superconducting 
temperature. The ratio Ar/ (ksT c ji) (with T c # the R— dependent transition temperature) is not constant any more 



4 



but depends on the width. For example, for the resonances at R — 0.86 nm and R = 1.12 nm (the panel (al) 
of FigQ we have found A, r /(1sbT c b) w 1.28. Hence, to get an idea about T c r/T c for 0.8 nm < R < 1.8 nm 
at n e = 3.878 nm~ 3 , one can simply multiply the data of FigQ] (al) by a factor of about 1.375. Situation with 
panel (bl) of Fig 2] is very similar. Here this factor is close to 1.66. Shape resonances in the dependence of the 
quasiparticle spectrum gap on the wire width are accompanied by strong increase and significant spatial variations 
of the order parameter A(p). Even for R > 5 nm the resonance effect is still considerable, and the superconducting 
order parameter differs significantly from the bulk value (see two examples of resonances in Fig|2(a)). In agreement 
with calculations for R > 10 nm Ll7j we observe strong oscillations of the order parameter with the frequency roughly 
proportional to kp. These oscillations are washed out at large radii about 100 nm (see the results of Ref. |l7|'). 

As it is shown, T c ^ sharply rises at a resonant point (more than by order of magnitude) and drops down to a 
value about T c between resonances. This makes it possible to expect that the effect will survive for fluctuating width 
in a smoothed variant depending on the scenario of fluctuations. Here it is necessary to mention the increase of the 
superconducting temperature by a factor of 1.1 for a single-crystal Sn nanowire with R = 10 nm Numerical 
investigation of Eqs. © and (Q) for the Sn parameters (gN(0) = 0.25 and hwo/ks = 195°K, see Ref. 0) indeed 
shows the presence of resonances with An/Abuik ~ 1-1 near the point R — 10 nm. One of them is situated at 
R = 9.642 nm, where A^/A^fe = 1.13. The superconducting order parameter at this resonant point has still 
significant spatial variations in the transverse direction and goes well above the bulk value (see Fig[2](b)). It is worth 
noting that when exploring the BdG equations for nanowires with R > 10 nm, the authors of Ref. |17| have found 
10% — 20% deviations of the averaged order parameter from the bulk limit even with radius fluctuations taken into 
account. 

In summary we investigated influence of the confinement on the transverse superconducting modes in a cylindrical 
nanowire at R < 10 nm in the clean limit. Numerical investigations of the BdG equations revealed a sequence of 
significant shape resonances in the dependence of the energy gap parameter on the nanowire width. The resonant 
deviations of An from the bulk value are by about order of magnitude larger than the shape resonances in the energy 
gap parameter of ultrathin metallic films of the same width. This makes it possible to expect that some smoothed but 
still profound effect will remain even in the presence of fluctuations of the wire radius (much less significant resonances 
in the superconducting order parameter survive for a nanowire with the average radius R ~ 10 nm for the variations 
SR < 0.5 nm [13). 

Acknowledgments 

We would like to thank A. Vagov for useful discussions and remarks. 



[1] L. P. Gor'kov, Sov. Phys. JETP 7, 505 (1958). 
[2] N. N. Bogoliubov, Sov. Phys. Usp. 67, 549 (1959). 
[3] L. N. Cooper, Phys. Rev. 104, 1189 (1956). 

[4] J. M. Blatt and C. J. Thompson, Phys. Rev. Lett. 10, 332 (1963). 

[5] M. Farine, F. W. Hekking, P. Schuck, and X. Vinas, Phys. Rev. B 68, 024507 (2003). 

[6] V. A. Ivanov, V. R. Misko, V. M. Fomin, J. T. Devreese, Solid State Commun. 125, 439 (2003). 

[7] W. Satula et al., Phys. Rev. Lett. 81, 3599 (1998). 

[8] S. Hilaire, J.-F. Berger, M. Girod, W. Satula, P. Schuck, Phys. Lett. B 531, 61 (2002). 

[9] M. Tian, J. Wang, J. S. Kurtz, Y. Liu, M. H. W. Chan, T. S. Mayer, T. E. Mallouk, Phys. Rev. B 71, 104521 (2005). 
[10] N. Giordano, Physica (Amsterdam) B 203, 460 (1994). 

[11] F. Sharifi, A. V. Herzog, and R. C. Dynes, Phys. Rev. Lett. 71, 428 (1993). 
[12] X. S. Ling et al, Phys. Rev. Lett. 74, 805 (1995). 

[13] A. D. Zaikin, D. S. Golubev, A. van Otterlo, G. T. Zimanyi, Phys. Rev. Lett. 78, 1552 (1997). 
[14] E. H. Hwang, S. Das Sarma, M. A. Stroscio, Phys. Rev. B 61, 8659 (2000). 
[15] P. G. de Gennes, Superconductivity of Metals and Alloys (W. A. Benjamin, New York, 1966). 
[16] A. L. Fetter, J. D. Walecka, Quantum Theory of Many-Particle Systems (Dover, Mineola, 2003). 
[17] J. E. Han, V. H. Crespi, Phys. Rev. B 69, 214526 (2004). 



